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I.  INTRODUCTION 


Military  specification  armor  steel,  herein  referred  to  as  rolled  homoge¬ 
neous  armor  (RHA),  is  important  in  ballistics  not  only  because  of  its  use  on 
combat  vehicles  to  resist  penetration  and  blast  damage,  but  also  as  a  target 
material  in  the  research  laboratory  to  evaluate  the  performance  of  new  penetra 
tor  designs  and  materials.  An  effective  evaluation  requires  a  knowledge  of 
penetrator-target  interactions,  which  include  the  initiation,  propagation  and 
interaction  of  large  amplitude  stress  waves  in  both  target  and  penetrator.  A 
knowledge  of  the  stresses,  strains,  particle  velocities  and  displacements  in 
the  target  material  would  aid  in  the  understanding  and  modeling  of  the  penetra 
tion  process.  This  understanding  then  could  be  applied  to  designing  the  con¬ 
figuration  and  placement  of  armor  on  combat  vehicles  and  fortifications. 

The  dynamic  mechanical  parameters  mentioned  above  also  should  be  related 
in  a  general  analytical  form,  that  is,  a  constitutive  equation,  as  this  form 
lends  itself  readily  and  with  flexibility  for  use  in  situations  involving 
combined  loading.  Thus  far,  the  only  analytic  relationships  among  the  mechani 
cal  properties  variables  for  RHA  are  for  elastic  deformations,1  2  where  Hooke' 
law  applies,  and  for  shock  loading  where  a  Hugoniot  equation  is  applicable. 3 

The  latter  was  deduced  from  plate  impact  experiments,  and  the  former  by 
an  ultrasonic  resonance  technique  and  tensile  stress-strain  tests  performed  at 
two  different  strain  rates.  Both  strain  rates  (e),  however,  were  such  that 
wave  propagation  effects  could  be  safely  ignored  (e  <  1  s-*) ,  and  little 
difference  in  material  properties  was  discernable  between  data  obtained  at  the 
different  rates.  Although  stress  and  strain  are  recorded  beyond  the  ultimate 
stress  at  both  strain  rates,  no  attempt  is  made  to  develop  an  analytical  model 
for  the  post-yield  deformation  behavior. 

The  objectives  of  this  study  of  the  deformation  of  RHA  are: 

•  obtain  compatible  uniaxial  stress  quasi-static  and  dynamic  post¬ 
yield  deformation  data 

•  examine  the  dynamic  data  as  a  wave  propagation  phenomenon,  and 
examine  both  types  of  data  in  the  context  of  extant  finite 
deformation  theory 


1 R .  F.  Bench,  "Quasi-Static  Tensile  Stress  Strain  Curves — II,  Rolled  Homoge¬ 
neous  Armor, "  BRL  Memorandum  Report  No.  2703,  Ballistic  Research  Laboratory , 
Aberdeen  Proving  Ground,  Maryland,  November  1976,  AD  #B016015L. 

2i?.  F.  Benck,  J.  L.  Robitaille,  "Tensile  Stress  Strain  Curves— III,  Rolled 
Homogeneous  Armor  at  a  Strain  Rate  of  0.42  s~^,"  BRL  Memorandum  Report  No. 
2760,  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland, 

June  1977. 

3G.  E.  Hauver,  "The  Alpha  Phase  Hugoniot  of  Rolled  Homogeneous  Armor, "  BRL 
Memorandum  Report  No.  2651,  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground,  Maryland,  August  1976,  AD  #B012871L. 
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•  correlate  data  obtained  over  a  wide  strain  rate  range  in  a  unified 
theory  of  post-yield  deformation  for  the  conditions  of  uniaxial 
stress . 

For  a  comparison  of  quasi-static  and  dynamic  data  to  be  meaningful,  the  deforma¬ 
tion  mode  (tension  or  compression)  must  be  identical  over  the  range  of  strain 
rates  examined.  Although  tensile  stress-strain  data  at  two  essentially  quasi¬ 
static  strain  rates  exists,  obtaining  tensile  wave  propagation  data  to  large 
strain  is  especially  difficult  for  high  strength  materials  such  as  armor  steel. 
Compression  experiments  thus  were  chosen. 

Quasi-static  stress-strain  curves  were  obtained  by  loading  resistance 
strain-gage  instrumented  samples  of  RHA  with  an  Instron  testing  machine 
equipped  with  a  data  logger,  which  recorded  both  stress  and  strain.  Dynamic 
loading  was  attained  by  subjecting  specimen  rods  to  axial,  symmetric,  free 
flight  impact  over  a  large  range  of  impact  velocities.  The  struck  rod  was 
instrumented  with  resistance  strain  gages  at  several  positions  along  its 
length,  and  the  gage  response  used  to  create  strain-time  profiles.  Additional 
data  obtained  from  the  impact  tests  were  duration  of  impact,  post  impact 
velocity  of  both  struck  and  striking  specimen,  and  residual  strain  profiles  of 
both  specimens. 

The  quasi-static  data  exhibited  Hookian  characteristics  up  to  a  well 
defined  yield  point,  and  post-yield  deformation  obeyed  Bell's  general  theory 
of  parabolic  plasticity.1*  The  parabola's  origin  was  at  the  yield  point,  and 
a  discrete  change  in  deformation  modulus  occurred  at  the  first  transition 
strain,  also  predicted  by  the  theory.  Dynamic  yield  stress  was  determined  by 
the  duration  of  impact  measurements,  a  sensitive  measure  of  the  dynamic  elastic 
limit  because  duration  of  impact  is  a  minimum  at  the  yield  velocity. 1+5  5  The 
dynamic  yield  stress  was  markedly  higher  than  the  quasi- static  value,  and  equal 
to  the  stress  at  the  first  transition  strain  of  the  quasi-static  stress-strain 
curve,  suggesting  a  heretofore  unexamined  relationship  between  dynamic  yield 
stress  and  the  transition  strain  structure  of  the  Bell  theory. 

Elastic  stresses  did  not  propagate  with  only  the  elastic  bar  wave  velocity, 
but  with  a  two-wave  structure  near  the  impact  end,  which  degenerated  to  the 
bar  velocity  at  distances  far  from  the  impact  end.  The  two-wave  structure 
corresponded  to  the  dilatational  and  shear  waves.  Levels  of  constant  strain 
beyond  yield  propagated  with  constant  velocity,  thus  establishing  the  applica¬ 
bility  of  the  one-dimensional  wave  propagation  theory  of  von  Karman,* * * 5 6  Taylor,7 


J.  F.  Bel l,  "A  Physical  Basis  fov  Continuum  Theories  of  Finite  Strain  Plastic¬ 

ity:  Part  I,"  Arch.  Rational  Mech.  Anal.,  70_,  1979 ,  pp  319-338,  and 

" .  Part  II,"  Arch.  Rational  Mech.  Anal.,  75,  1981,  pp  104-126. 

5J.  F.  Bell,  "On  the  Dynamic  Elastic  Limit, "  Exp.  Mech.,  22,  No.  7,  1982, 
pp  270-276. 

6T.  von  Karman,  "On  the  Propagation  of  Plastic  Deformation  in  Solids, "  National 
Defense  Research  Council  Report  A-29,  OSRD  365,  U.S.A.,  February  1942. 

7G.  I.  Taylor,  "The  Plastic  Wave  in  a  Wire  Extended  by  an  Impact  Load,"  Civil 
Defense  Research  Committee  Report  R.C.  329,  British  Ministry  of  Rome  Security, 
1942.  8 


Rakhmatulin,& * 8  and  White  and  Griffis.9  The  plastic  wave  speeds,  Cp  obtained 

were  given  by  Cp  =  /(9a/9e) /p  with  9a/9e  the  slope  of  the  static  stress-strain 

curve,  and  p  the  mass  density.  The  one  dimensional  wave  propagation  theory 
also  predicts  an  integral  relationship  between  particle  velocity  u  and  strain 

e:  u  =/e  C  de,  where  C  is  the  plastic  wave  speed.  The  symmetry  of  the  impact 
P  P  1 

allows  an  accurate  determination  of  maximum  particle  velocity,  and  the  measured 
strain  maxima  near  the  impact  end  agree  closely  with  prediction  when  the  wave 
speed  expression  as  given  above  is  the  integrand. 


The  post-yield  deformation  behavior  of  rolled  homogeneous  armor  has  been 
characterized  by  a  parabolic  constitutive  equation  applicable  over  a  strain 
rate  range  of  several  orders  of  magnitude.  The  quasi-static  and  dynamic  stress- 
strain  curves  are  not  identical  but  are  related  through  a  general  theory  of 
plasticity. 


II.  THEORETICAL  BACKGROUND 

A.  Parabolic  Response  Function  for  Plastic  Deformation 

Data  obtained  from  over  2,000  large  deformation  tests  on  27  annealed 
metals  and  16  "as  received"  structural  metal  alloys  of  aluminum,  copper 
and  iron  exhibit  plastic  constitutive  equations  whose  general  form  is 
parabolic.4*10’11’12’13  The  full  details  of  the  theory  and  supporting  experi¬ 
mental  evidence  can  be  found  in  the  literature,  but,  for  ready  reference,  the 
most  important  aspects  pertaining  to  this  study  of  uniaxial  stress  loading  are 
included  here. 


&K.  A.  Rakhmatulin ,  "Propagation  of  a  Wave  of  Unloading, "  Soviet  J.  Appl. 

Math.  Mech. ,  9_,  1945. 

9 M .  P.  White,  Le  van  Griffis,  "The  Permanent  Strain  in  a  Uniform  Bar  Due  to 

Longitudinal  Impact,"  J.  Appl.  Mech. ,  14,  1947,  pp  A-537  -  A-342. 

10J.  F.  Bell,  The  Physics  of  Large  Deformation  of  CrpystaJ/line  Solids ,  Springer 
Tracts  in  natural  Philosophy,  Vol.  14,  Spring er-Ver lag ,  Heidelberg,  Berlin, 
New  York,  1968. 

Dj.  f.  Bell,  The  Exjpervnerital  Foundations  o$_  Solid  Mechanic^,  Handbuch  der 
Physik ,  Vol.  VIa/1,  Spring er-Ver lag,  Heidelberg,  Berlin,  New  York,  1973. 

12 j.  F.  Bell,  "Technological  Perspectives  from  Two  Decades  of  Fundamental 
Research  in  Dynamic  Plasticity,"  BRL  Contract  Report  No.  184,  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland ,  October  1974, 

AD  M003175. 

13J.  F.  Bell,  " Origins  in  Experiment  of  a  New  General  Theory  of  Plasticity  for 
Structural  Metal  Alloys,  "  BRL  Contract  Report  No.  311,  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland,  August  1976,  AD  #A029845. 
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(1) 


Early  studies  by  Bell  on  annealed  aluminum  show  a  plastic  engineering 
stress -engineering  strain  relationship  given  by 


where  g  -  parabola  coefficient,  GPa 
Cp  -  engineering  strain,  mm/mm 
Op  -  engineering  stress,  GPa 

and  subscript  p  refers  to  plastic  conditions.  The  parabola  coefficient  varies 
linearly  with  temperature  and  depends  on  the  value  of  the  material's  isotropic 
linear  elastic  shear  modulus  at  absolute  zero.  It  is  given  by 

3  =  (2/3)r/2  pCO)  Bq  (1-Ta/Tm)  (2) 

where  r  -  non-negative  integer  mode  index 

Bq  -  dimensionless  universal  constant,  0.0280 

T&  -  material  temperature,  degrees  Kelvin 

Tm  -  material  melting  temperature,  degrees  Kelvin 

y(0)  -  isotropic  linear  elastic  shear  modulus  at  zero  degrees 
Kelvin,  GPa 

2 

The  parabolicity  is  obvious  when  the  stress-strain  data  are  viewed  as  a  vs *  *  e 
plots,  as  the  slopes  become  linear.  A  similar  parabolic  response  function 
applies  to  other  metals.14 

2 

The  slope  of  a  o  vs.  e  plot  need  not  be  single  valued,  but  may  have 
discrete  values  corresponding  to  discrete  changes  in  the  mode  index  r,  which 
may  either  increase  or  decrease.  The  changes  in  slope  occur  at  specific 
values  of  plastic  strain  given  by 

eN  -  (2/3)N/2//3  (3) 

where  N  may  have  the  values  18,  13,  10,  8,  6,  4,  2,  and  0,  in  the  order  of 
increasing  plastic  strain.  These  eN  are  transition  strains,  and  although  the 

physics  is  not  yet  explicable,  several  other  physical  properties  exhibit 
radical  changes  when  materials  are  deformed  beyond  those  values.11’* 


14J.  F.  Belt ,  "Generalized  Large  Deformation  Behaviour  for  Faoe-Centered- Cubic 
Solids:  Bickel ,  Aluminum ,  Gold ,  Silver  and  Lead,"  Phil.  Mag.  11,  No.  114, 
1965,  pp  1135-1156. 

*J •  F.  Belt,  "The  Initial  Yield  and  Large  Deformation  of  Mild  Steel  in  the 
Context  of  a  General,  Incremental  Theory  of  Finite  Plastic  Strain"  (manu¬ 
script  pending  publication) . 


Bell10  shows  an  empirical  temperature  dependent  relationship  for  the 
linear  elastic  shear  modulus 

p(Ta/Tm)  =  1.03  (l-Ta/2Tm)  y(0)  (4) 

which  is  valid  over  the  fractional  temperature  range  0.06  £  T  /T  i  300 /T. 

ci  III  III 

This  observation,  along  with  the  elastic  constant  relation 

v  =  E/2 (1+v)  (5) 


where  E  -  Young's  modulus,  GPa 
y  -  shear  modulus,  GPa 
v  -  Poisson's  ratio,  dimensionless 
was  used  to  modify  Eq.  (2)  in  terms  of  E  and  v: 

3  =  (2/3)r/2  Bq  {[E/1.03(l+v)][l-(2-Ta/Tm)'1]}  (6) 


The  bracketed  portion  of  Eq.  (6),  designated  by  { }  thus  becomes  the  material 
property  constant  in  the  parabola  coefficient. 

For  fully  annealed  materials  with  a  negligible  elastic  limit,  the  total 
stress  and  total  strain  may  be  used  in  Eq.  (1).  Materials  with  an  appreciable 
elastic  limit  also  exhibit  parabolic  plasticity,  but  the  origin  of  the  parabola 
must  be  shifted  to  account  for  the  elastic  region.  The  counterpart  of  Eq.  (1) 
is  thus 

a  -  oy  =  Bfc-Ey)'2  (7) 


where  a  and  e  refer  to  total  nominal  stress  and  total  nominal  strain  respec¬ 
tively,  and  the  subscript  y  refers  to  values  at  yield.  This  region  of  deforma¬ 
tion  where  reference  is  made  to  the  elastic  region  is  termed  the  intermediate, 
or  elastic-plastic  region.  The  parabola  coefficient  B  has  the  same  initial 
mode  index  as  the  annealed  material,  and  changes  in  mode  index  occur  in  the 
intermediate  region  as  well.  The  displacement  of  the  parabola's  origin  in 
stress-strain  space,  however,  must  be  accounted  for  in  computing  the  transition 
strains,  thus 


where  e"  is  the  total  strain  at  which  a  transition  can  occur. 

N 

B .  Dynamic  Yield  Condition 

The  yield  stress  and  strain  for  quasi-static  loading  can  be  determined  by 
locating  the  terminus  of  the  linear  elastic  region  of  the  stress-strain  curve, 
or  in  cases  where  no  definitive  yield  point  occurs,  by  the  0.2  percent  strain 
offset  method.  In  either  case,  load  and  deformation  can  be  measured  directly 
and  related  to  stress  and  strain.  For  dynamic  loading,  one  may  use  the  linear 
elastic  energy  balance  equation 
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(9) 


1/2  ae  =  1/2  pif 

and  the  stress-strain  relation  for  uniaxial  loading 

a  =  XA  +  2ye  (10) 

[or  a  =  Ee  (1C)')] 

where  u  -  particle  velocity,  mm/ys 

A  -  Lame’s  constant,  GPa 

3 

P  -  mass  density,  g/cm 
3  3 

A  -  dilatation,  mm  /mm 

to  obtain  a  relation  between  stress  and  particle  velocity 

a  =  pCQu  (11) 

where  -  i/e/ p ,  the  elastic  bar  wave  speed,  mm/ys.  This  relation  is  valid 

for  u  <  u  ,  that  is  for  particle  velocities  below  the  elastic  limit  velocity, 
and  hence  for  elastic  stresses. 

Bell  *  shows  that  a  measure  of  the  duration  of  impact  for  axial,  free 
flight  collisions  of  two  identical  long  rod  specimens  is  a  sensitive  indicator 
of  the  end  of  the  elastic  region.  Sears  shows  that  duration  of  impact 
decreases  as  impact  velocity  increases  in  the  elastic  region.15  Bell  shows 
the  converse  is  true  in  the  plastic  region,  and  that  a  minimum  occurs  at  the 
transition  from  elastic  to  plastic  behavior.5  The  maximum  particle  velocity 
at  the  minimum  duration  of  impact  is  u  ,  which  can  be  used  in  Eq.  (11)  to 
calculate  dynamic  yield  stress.  ^ 

Another  determination  of  the  particle  velocity  at  the  end  of  the  elastic 
region  can  be  made  by  examining  the  coefficient  of  restitution  versus  impact 
velocity.  The  coefficient  e,  is  defined  as  the  absolute  value  of  the  ratio  of 
the  final  relative  velocity  of  the  specimen  rods  to  their  initial  relative 
velocity.  It  remains  essentially  constant  for  elastic  deformations,  but  begins 
a  sharp  decrease  when  some  of  the  impact  energy  is  used  to  produce  plastic 
deformation,  that  is,  at  yield. 

C .  Plastic  Wave  Propagation 

Plastic  deformation  of  solids  by  stress  waves  is  a  complex  process. 
Strain-time  measurements  on  long  rod  specimens  subjected  to  axial,  symmetric, 
constant  velocity  impact  provide  experimental  data  amenable  to  a  one-dimensional 
analysis .  The  analysis  used  is  the  Lagrangian  form  of  a  strain-rate  indepen¬ 
dent,  finite  amplitude  wave  propagation  theory. ^  The  development  of  the  theory 


15J.  E.  Sears,  "On  the  Longitudinal  Impact  of  Metal  Rods  with  Rounded  Ends, 
Part  II,"  Cambridge  Phil.  Soc.,  21_,  1912,  pp  49-105. 
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from  first  principles  is  reproduced  by  Kolsky, so  only  the  most  important 
features  are  shown  here.  The  salient  feature  is  that  each  level  of  strain 
propagates  with  a  velocity  dependent  upon  the  instantaneous  slope  of  the 
governing  stress-strain  curve,  which  need  not  be  specified  beforehand.  This 
feature  arises  from  the  solution  of  a  wave  equation 


P 


(12) 


where  t  -  time,  ys 

u  -  displacement,  mm 

S  -  deformation  modulus,  3 a/3e,  GPa 

X  -  Lagrangian  coordinate,  mm. 

(This  notation  parallels  that  used  by  Kolsky.16) 

The  deformation  modulus  S  need  not  be  initially  specified,  but  must  be  a  single 
valued  function  of  strain,  S(e) ,  and  the  stress-strain  curve  must  be  concave 
toward  the  strain  abscissa.  The  stress  wave  propagation  velocity  for  levels 
of  constant  strain  is  given  by 


C(e)  =  /S(E)/p  (13) 

That  levels  of  constant  strain  propagate  with  constant  velocity  may  be  experi¬ 
mentally  determined  by  comparing  wave  speeds  for  levels  of  strain  over  at  least 
two  portions  of  a  rod  subjected  to  a  constant  velocity  impact.  This  is  the 
first  requirement  to  determine  the  applicability  of  the  theory. 

The  second  requirement  of  the  theory  is  that  a  relationship  between 
particle  velocity  and  corresponding  strain  exist.  Integrating  Eq.  (12) 
produces 


u  =  /£C(e)  de  (14) 

which  may  also  be  checked  experimentally  by  comparing  maximum  strains  and 
maximum  particle  velocities.  It  follows  that  the  stress  is  given  by 

a  =  pJeC^(e)  de  (15) 

The  wave  equation  [Eq.  (12)]  may  be  integrated  graphically,  numerically,  or  if 
the  wave  speeds  are  experimentally  found  to  be  functionally  related  to  the 
strains,  analytically.  The  results  are  experimentally  determined  particle 
velocity-strain  and  dynamic  stress-strain  curves. 


1  Kolsky ,  Stress  Waves  in  Solids,  Dover  Publications ,  Inc.,  New  York,  1963, 

pp  164-170. 


13 


III.  EXPERIMENTAL 


A.  Specimen  Material  and  Preparation 

Rolled  homogeneous  armor  is  a  high  nickel,  chrome-moly  steel.  Chemical 
analysis  of  the  material  used  is  presented  in  Table  1. 


TABLE  1 

CHEMICAL  ANALYSIS  OF3  38  mm  RHA 


Element 

Weight  Percent 

Iron 

>  94.82 

Nickel 

3.04 

Chromium 

1.07 

Molybdenum 

0.10/0.25 

Manganese 

0.27 

Carbon 

0.26 

Silicon 

0.18 

Copper 

0.07 

Aluminum 

<0.02 

Vanadium 

<  0.01 

Sulphur 

0.008 

Phosphorous 

0.001 

100.00 


After  having  been  cast  into  ingots,  the  steel  is  hot  rolled  and  cross-rolled 
to  various  standard  plate  thicknesses,  and  heat  treated,  quenched  and  tempered 
in  accordance  with  MIL  SPEC  MIL  12560B.17  The  specification  requires  that 
hardness  be  kept  within  certain  limits  depending  on  plate  thickness.  The 
hardness  of  the  thirty-eight  millimeter  thick  plate  stock  from  which  specimens 
were  prepared  was  30-31  Rockwell  C,2  although  a  cyclic  hardness  variation 
across  the  thickness  of  the  plate  has  been  reported.3 

The  plate  stock  was  saw  cut  into  bars  38  x  38  x  300  mm.  The  bar's 
longitudinal  axis  was  in  the  direction  of  the  final  roll  of  the  plate,  which 
was  determined  by  a  macroscopic  examination  of  the  flow  lines  after  etching 
with  Frye's  reagent.18  The  choice  of  axis  orientation  relative  to  the  plane 


1 7 Military  Specification ,  Steel  Armor:  Plate  Wrought  Homogeneous;  Combat 
Vehicle  Type  (1/4  to  6  inches ,  vncU) ,  MIL  SPEC  MIL-S-12560B  (ORD) ,  31  July 
1962. 

18 G.  L.  Kehl ,  Principles  of  Metallo^raphic  Laboratory  Practice ,  McGraw-Hill 
Book  Co.,  Inc.,  New  York,  1949,  pp  205-206,  Table  28. 
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of  the  plate  was  arbitrary,  but  a  single  orientation  of  all  test  specimens 
prevented  variations  in  the  data  due  to  preferred  orientation  of  the  grain 
structure  induced  by  rolling.  The  bars  were  machined  to  flat- ended  cylindrical 
rods  with  a  surface  finish  of  sixty-three  microinch. 

The  quasi-static  compression  specimens  were  9.55  mm  in  diameter  and 
28.58  mm  long,  providing  a  length  to  diameter  ratio,  (L/D)  of  three  to  one. 

The  foil  type  resistance  strain  gages  used  to  measure  strain  were  high  elonga¬ 
tion,  90  degree  Tee  rosettes,  with  1.57  mm  gage  length.  Each  specimen  was 
instrumented  with  three  gages,  120  degrees  apart,  midway  along  the  length  of 
the  specimen.  The  two  gage  axes  of  the  rosette  were  parallel  and  perpendicu¬ 
lar,  respectively,  to  the  specimen  axis. 

The  long  rod  impact  specimens  were  25.10  mm  diameter  and  249.75  ±  4.45  mm 
long,  for  an  L/D  of  ten.  The  impact  end  of  each  rod  was  ground  in  a  jig  to 
insure  its  perpendicularity  to  the  longitudinal  axis  of  the  rod  and  to  provide 
a  flat  surface  to  facilitate  axial  alignment  of  the  striker  and  struck  rods 
prior  to  impact.  Rods  were  paired  for  each  test  using  the  criteria  that  the 
difference  in  their  lengths  be  0.03  mm  or  less.  The  struck  specimens  were 
instrumented  with  foil  type  resistance  strain  gages.  A  pair  of  3.18  mm  gage 
length  gages,  180  degrees  apart,  were  applied  25.0,  50.0,  and  75.0  mm  from  the 
impact  end,  with  the  gage  axis  parallel  to  that  of  the  rod.  A  pair  of  parallel, 
circumferential  scribe  lines,  approximately  10  mm  apart,  straddled  each  strain 
gage  position.  The  pre-  and  post- impact  line  spacings  were  used  to  determine 
residual  strains.  The  striker  specimens  also  had  scribe  line  pairs  at  the 
same  positions  as  the  struck  specimens.  The  scribe  line  and  gage  center 
position  measurements  were  made  using  a  travelling  stage  microscope. 

B .  Experimental  Procedures 

The  quasi-static  compression  tests  were  performed  with  a  standard  Instron 

testing  machine  at  a  strain  rate  of  3.0  x  10  ^  s  The  instantaneous  resist¬ 

ance  R,  of  each  resistance  strain  gage  was  monitored  and  strain  calculated  by 

e  =  /r7R0  -  1  (16) 

where  R  is  the  initial  resistance  of  the  unstrained  strain  gage.  This 
o 

expression  for  strain  results  when  a  gage  factor  of  (2  +  e)  is  used.  The 
details  of  the  apparatus  and  data  reduction  can  be  found  elsewhere, *9  however, 
measured  strains,  both  longitudinal  and  lateral,  were  accurate  to  ±  0.0006  per¬ 
cent  strain. 

The  experimental  apparatus  used  for  the  constant  velocity  impact  tests  is 
schematically  illustrated  in  Figure  1.  The  smooth  bore  gas  gun  had  been  honed 
to  an  i.d.  of  25.15  mm.  Compressed  dry  nitrogen  gas  was  the  propellant,  and 
firing  actuated  by  a  solenoid  valve.  The  muzzle  of  the  1.83  m  long  barrel 
was  sufficiently  ported  with  several  10  mm  diameter  holes  to  relieve  the 


l3R.  F.  Bench,  D.  A.  DiBerardo,  R.  E.  Franz ,  "Quasi- Static  Compression  Tests  - 
S7  Toot  Steel, rf  BRL  Memorandum  Report  No,  0367 ,  Ballistic  Research  Laboratory , 
Aberdeen  Proving  Ground,  Maryland,  October  1980 . 
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propellant  gas  pressure  so  that  at  impact  the  striker  was  in  essentially  free 
flight  with  constant  velocity.  The  impact  chamber,  which  enclosed  the  muzzle 
and  struck  specimen,  was  evacuated  to  less  than  100  yin  of  mercury  pressure. 

The  vacuum  eliminated  an  air  "cushion"  between  the  flat  impact  ends.  The  exit 
port  of  the  chamber,  which  led  to  a  catch  tank  cushioned  with  textile  waste, 
was  sealed  with  a  mylar  diaphragm. 

The  struck  specimen  was  supported  in  a  cantilever  arrangement  with  four 
#10-32  nylon  screws.  The  hemispherical  tip  screws  were  threaded  into  the  speci¬ 
men  holder,  and  contacted  the  specimen  far  from  the  impact  end.  Axial  align¬ 
ment  was  achieved  by  mating  the  ground  impact  ends  of  the  striker  and  struck 
specimen  in  the  position  they  would  occupy  at  impact.  (Approximately  one-half 
of  the  striker's  length  was  in  the  muzzle).  The  struck  specimen's  angular 
position  was  adjusted  with  the  nylon  screws  until  no  light  could  be  observed 
between  the  specimens,  indicating  total  areal  contact  of  the  flat  impact  ends. 
The  striker  was  then  drawn  to  the  breech  of  the  gas  gun. 

The  optical  system  used  to  measure  striking  velocity,  duration  of  speci¬ 
men  contact,  and  post- impact  velocity  of  struck  and  striker  specimens  was 
analogous  to  that  which  Bell21"1  uses.  The  essential  components  were  a  collimated 
light  field,  which  illuminated  a  velocity  field,  and  a  photomultiplier  tube 
(PMT) •  A  small  filament  DC  incandescent  lamp,  mounted  at  the  focus  of  a  lens 
comprised  the  light  system.  Light  field  size  was  controlled  by  opaque  baffles. 
The  velocity  field,  which  was  an  integral  part  of  the  specimen  holder,  con¬ 
sisted  of  a  narrow  open  slot,  6.5  x  100  mm,  across  which  two  pairs  of  2.5  mm 
wide  metal  strips  had  been  fastened  with  screws.  The  plane  of  the  strips  was 
approximately  8  mm  below  the  lowest  portion  of  the  specimens.  The  illuminated 
area  was  painted  flat  black  to  minimize  reflections  which  could  cause  erroneous 
signals.  Light  that  passed  through  the  field  slot  was  collected  by  a  con¬ 
densing  lens  and  impinged  on  the  collector  of  the  PMT .  Output  voltage,  which 
depended  on  the  amount  of  light  reaching  the  collector,  was  monitored  with  an 
oscilloscope. 

As  the  striker  approached  the  struck  specimen,  it  covered  the  velocity 
field  and  less  light  reached  the  PMT,  causing  a  decrease  in  output  voltage. 

The  light  level  remained  constant  as  an  opaque  strip  was  traversed,  as  did 
output  voltage.  This  constant  output  was  easily  detected  as  a  "flat"  on  the 
oscillogram,  typified  by  Figure  2.  The  time  to  traverse  the  known  distance 
between  the  strips  was  measured  with  the  aid  of  the  50  ys  timing  signal,  and 
striking  velocity  was  determined.  When  the  specimens  collided,  the  velocity 
field  was  completely  covered  and  no  light  reached  the  PMT.  The  "dark"  condi¬ 
tion  remained  as  long  as  the  specimens  were  in  contact.  As  the  specimens 
rebounded  and  separated,  the  struck  specimen  passed  over  a  second  pair  of 
opaque  strips.  This  action  produced  the  rising  signal  portion  of  the  Type  I 
trace,  and  the  post-impact  velocity  of  the  struck  specimen  was  determined. 

The  striker  specimen  too,  passed  over  the  second  pair  of  strips  a  short  time 
later,  as  can  be  seen  on  the  upper  trace  of  Type  II,  where  the  recording  time 
was  ten  times  that  of  the  Type  I  signal.  The  lower  trace  in  Type  II  is  that 
portion  of  the  PMT  signal  during  specimen  contact,  recorded  on  a  faster  time- 
base  to  improve  resolution.  Thus,  time  measurements  on  the  PMT  signal  provided 


20J.  F.  Bell ,  "An  Experimental  Study  of  the  Unloading  Phenomenon  in  Constant 
Velocity  Impact J.  Mech.  Phys.  Sol.,  9j  1961,  pp  1-15. 
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Figure  1.  Schematic  diagram  of  impact  test  apparatus;  the  optical 
velocity  measuring  field  is  shown  in  the  inset. 
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Figure  2.  Typical  PMT  signals.  Type  I  was  recorded  on  time  base  I.  Type  II 
was  delayed  and  recorded  on  time  bases  ten  times  longer  (upper) , 
and  one-tenth  as  long  (lower),  as  time  base  I. 


all  information  needed  to  compute  duration  of  impact  and  coefficient  of 
restitution. 

Time  resolution  was  limited  by  the  resolution  of  the  digitizer  used  for 
distance  measurements  on  the  oscillograms,  and  the  relatively  long  times  for 
specimens  to  traverse  the  velocity  field.  Time  resolution  on  the  PMT  oscillo¬ 
grams  was  typically  ±  4.1  microseconds,  and  was  the  limiting  factor  in  deter¬ 
mining  velocity,  duration  of  impact,  and  coefficient  of  restitution.  Striking 
velocities  were  accurate  to  within  ±  0.68  meters  per  second.  Duration  of 
impact  and  coefficient  of  restitution  values  were  accurate  to  within  ±  3.2  per¬ 
cent  and  ±  2.4  percent  respectively,  of  their  determined  values. 

Dynamic  strain  measurements  were  made  with  foil  type  resistance  strain 
gages,  positioned  as  described  above.  Each  gage  was  powered  by  an  individual 
capacitive  discharge  power  supply  through  a  bridge  resistance  network.  The 
current  was  pulsed  on  several  microseconds  before  impact  and  shunted  away  from 
the  strain  gage  after  the  recording  time  interval,  typically  100  microseconds. 
Output  voltage  was  monitored  on  an  oscilloscope  and  recorded  on  photographic 
film.  Calibration  of  strain  gages  and  recording  devices  was  performed  before 
each  test  by  introducing  resistances  in  parallel  with  the  strain  gage  and 
pulsing  the  power  supply.  The  calibration  resistances  used  simulated  com¬ 
pressive  strains  of  approximately  one  and  five  percent.  The  output  voltages 
obtained  during  calibration  were  used  to  calculate  a  voltage  deflection  factor. 
The  instantaneous  resistance  of  the  strain  gage  was  determined  by  solution  of 
the  applicable  bridge  network  equations,  and  strain  was  computed  by  Eq.  (16) . 
Dynamic  strain  measurements  were  estimated  to  be  accurate  to  within  ±0.05 
percent  strain. 

Wave  speeds  for  levels  of  constant  strain  were  determined  from  the  char¬ 
acteristics  of  a  Lagrangian  diagram  constructed  by  plotting  arrival  time  of  a 
level  of  strain  versus  distance  from  the  impact  end.  An  external  estimate  of 
the  uncertainty  associated  with  the  wave  speed  determinations  was  made.  The 
uncertainty  was  approximately  ±  0.12  mm/ys.  Wave  speeds  found  to  be  constant 
within  this  limit  were  assumed  to  be  constant  throughout,  and  the  character¬ 
istic  was  determined  by  least  squares  methods. 

Residual  strain  at  positions  along  the  specimen  rods  was  determined  by 
two  methods.  Method  one  utilized  the  pre-  and  post- impact  spacing  of  the 
circumferential  scribe  lines  described  above.  The  residual  strain  so  deter¬ 
mined,  although  distributed  over  a  finite  gage  length,  was  associated  with  a 
rod  position  midway  between  the  two  scribe  lines  used  to  determine  it.  Method 
two  used  the  pre-  and  post-impact  diameters,  which  were  measured  with  a  blade 
micrometer,  to  determine  diametral  strain.  Diametral  strain  was  converted  to 
longitudinal  strain  by  assuming  isochoric  deformation.  The  residual  strains 
determined  by  both  methods,  for  both  the  striker  and  struck  specimens  were 
averaged  for  each  position.  These  values  were  estimated  to  be  accurate  to 
within  ±  0.024  percent  strain. 
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IV.  RESULTS 


A.  Quasi-Static  Compression  Tests 

The  stress-strain  curves  for  the  two  quasi-static  compression  tests  are 
presented  in  Figure  3.  The  initial  portion  of  the  loading  curve  was  linear 
with  the  proportional  limit  at  0.812  GPa  stress  and  0.391  percent  strain. 
Young's  modulus  is  208.0  GPa,  Poisson's  ratio  0.290,  and  the  shear  modulus  is 
computed  to  be  80.6  GPa. 

A  general  power  function  equation 

a  -  cy  =  H(e  -  ey)a  (17) 

was  assumed  for  the  post-yield  stress-strain  data  in  nominal  form,  where  H  is 
the  power  function  coefficient  and  a  is  a  dimensionless  exponent.  A  plot  of 
the  data  as  the  common  logarithm  of  each  side  of  Eq.  (17),  where  0.821  GPa  and 
0.404  percent  strain  were  used  as  the  yield  stress  and  strain  (Figure  4),  was 
linear  with  the  slope  (a)  equal  to  0.517,  indicating  a  nearly  parabolic  func¬ 
tion  . 

B.  Dynamic  Yield  Stress 

Fifteen  symmetric,  free  flight,  axial  impact  experiments  were  performed, 
with  striking  velocities  ranging  from  20  m/s  to  over  100  m/s.  Time  of  specimen 
contact  during  impact  and  coefficient  of  restitution  were  determined.  Figure  5 
shows  the  measured  time  of  contact  (divided  by  specimen  length  to  normalize 
variations  in  rod  lengths)  versus  maximum  particle  velocity,  which  is  one-half 
the  striking  velocity.  The  minimum  duration  of  impact  occurs  at  a  particle 
velocity  of  0.0236  ±  .0003  mm/ys.  The  coefficient  of  restitution  e,  versus 
maximum  particle  velocity  is  plotted  in  Figure  6.  It  is  essentially  constant 
(0.932  ±  .011)  up  to  the  particle  velocity  associated  with  the  duration  of 
impact  minimum,  after  which  it  decreases  non-1 inear ly  with  increasing  velocity. 
No  permanent  deformation  was  detected  on  specimen  rods  struck  below  47.2  m/s 
(u  =  0.0236  mm/ys). 

Utilizing  this  value  for  u,  along  with  the  quasi-static  value  of  Young's 
modulus  and  a  density**  of  7.848  g/cm^  in  Eq.  (11),  the  dynamic  yield  stress 
is  computed  to  be  0.954  *  .014  GPa.  The  dynamic  yield  strain,  computed  with 
Eq.  (107),  is  0.458  *  .008  percent.  This  value  may  be  compared  with  the 
magnitude  of  the  elastic  precursor  obtained  from  the  resistance  strain  gage 
measurements.  Figure  7  shows  strain-time  profiles  obtained  at  positions  50  mm 
and  ,75  mm  from  the  impact  end  of  the  rod,  at  three  impact  velocities.  In 
Table  2  are  the  maximum  particle  velocities  and  average  plastic  strain  rates 
for  the  data  shown  in  Figure  7. 


** Private  communication  from  R.  F.  Benck ,  Ballistic  Research  Laboratory . 
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Figure  3.  Stress  vs.  compressive  strain  at  a  strain  rate  of  3.0 
compared  with  parabolic  response  function. 
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Figure  4.  Log  post-yield  stress  vs.  log  post-yield  strain  for  quasi-static 
data.  Yield  stress  taken  as  0.821  GPa,  and  yield  strain  0.0040. 
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Figure  5.  Normalized  time  of  contact  vs.  maximum  particle  velocity. 
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Figure  6.  Coefficient  of  restitution  vs.  maximum  particle  velocity. 
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Figure  7.  Compressive  strain  vs.  time  after  impact  at  two  positions 
on  struck  rod  for  three  impact  velocities. 

25 


TABLE  2 


MAXIMUM  PARTICLE 

VELOCITY  AND  AVERAGE 

STRAIN  RATE 

Test  # 

“max 

S  (s'1’ 

17 

0.0319 

65.4 

20 

0.0249 

14.4 

22 

0.0515 

235.3 

Although  the  plastic  strain-rate  range  is  larger  than  an  order  of  magni¬ 
tude,  the  computed  value  of  yield  strain  is  in  excellent  agreement  with  the 
magnitude  of  the  elastic  precursor,  regardless  of  strain  rate. 

C .  Elastic  Wave  Propagation 

Examination  of  the  elastic  precursor  portion  of  the  strain  time  profiles 
of  Figure  7  showed  that  the  rise  to  yield  strain  was  not  linear.  Wave  speeds 
of  pre-yield  strains  were  determined,  and  two  features  were  evident:  (1)  a  two- 
wave  structure  was  present,  and  (2)  elastic  strains  did  not  propagate  with 
constant  velocity  along  the  rod.  Figure  8  shows  wave  speeds  measured  in  two 
regions  of  the  rod:  25-50  mm  and  50-75  mm  from  the  impact  end.  Also  indicated 
are  the  dilatational  wave  velocity,  given  by  /(A+2y)/p,  the  shear  wave 

velocity,  given  by  /y/p,  and  the  elastic  bar  wave  velocity  CQ.  The  first 

value  is  the  average  obtained  from  plate  impact  experiments, 3  5.83  mm/us,  and 
the  latter  two  computed  using  the  quasi-static  values  of  y  and  E,  and  are 
3.21  mm/ys  and  5.15  mm/ys,  respectively.  In  the  region  of  the  rod  nearer  the 
impact  end,  strains  below  approximately  0.20  percent  strain  propagated  with 
the  dilatational  velocity.  Strains  above  this  value,  up  to  the  dynamic  yield 
strain,  propagated  with  approximately  the  shear  wave  velocity.  The  data  of 
the  farther  region  of  the  rod  displayed  a  similar,  although  less  distinct  two- 
wave  structure.  The  change  in  wave  speeds  occurred  at  the  same  value  of  strain 
as  in  the  region  nearer  the  impact  end,  but  the  wave  speeds  were  more  distrib¬ 
uted.  The  statistical  mode  wave  speeds  for  the  strain  regions  above  and  below 
0.20  percent  are  =  5.58  mm/ys  and  C2  =  3.95  mm/ys. 

Lame's  constant  X  can  be  determined  from  knowledge  of  the  dilatational 
wave  speed  and  is  computed  to  be  105.5  GPa.  It  follows  that  the  bulk  modulus  K, 
is  159.3  GPa,  and  the  bulk  velocity  of  4.51  mm/ys  is  in  agreement  with  pre¬ 
viously  reported  values.3 

D .  Plastic  Wave  Propagation 

Levels  of  strain  greater  than  the  dynamic  yield  strain  propagated  with 
constant  velocity  beyond  the  first  diameter  length  from  the  impact  end. 
Increasing  levels  of  strain  propagated  with  slower  velocities.  The  constancy 
of  propagation  velocity  could  not  be  determined  for  levels  of  strain  which  did 
not  travel  to  all  three  strain-gage  positions,  so  plastic  wave  speed  data  were 
limited  to  strain  levels  of  approximately  1.4  percent  strain.  Single  wave 
speeds  for  levels  of  constant  strain,  and  decreasing  wave  speeds  for  increasing 
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Wave  speed  vs.  compressive  elastic  strain  for  two  regions  of  the  struck  rod. 


strains  indicated  that  the  first  applicability  requirement  of  the  strain-rate 
independent,  one  dimensional,  finite  amplitude  wave  propagation  theory  had  been 
met . 

A  general  power  function  equation  [Eq.  (17)]  was  assumed  as  the  governing 
stress-strain  relationship,  and  differentiated  to  obtain  an  expression  for  wave 
speeds : 


(18) 


where  e  is  the  dynamic  yield  strain.  The  numerator  under  the  square  root  sign 

corresponds  to  the  deformation  modulus  S(e),  in  Eqs.  (12)  and  (13)  .  A  loga¬ 
rithmic  plot  of  experimentally  determined  wave  speeds  versus  post-yield  strain 
is  shown  in  Figure  9.  The  relationship  is  linear  with  the  exponent  a  .=  0.522, 
indicating  a  nearly  parabolic  stress-strain  relationship,  as  was  found  for  the 
quasi-static  post-yield  deformation. 


V.  DISCUSSION  OF  RESULTS 


The  experimental  evidence  indicating  a  nearly  parabolic  stress -strain 
relation  suggests  examination  of  the  data  in  the  context  of  Bell's4  general 
theory  of  parabolic  plasticity.  The  material  has  an  appreciable  elastic 
limit,  therefore  the  elastic-plastic  equation  [Eq.  (7)3*  if  any>  should  apply. 
The  parabolic  response  function  was  postulated  and  made  linear  with  respect  to 
strain.  Quasi-static  yield  stress  was  determined  by  using  the  value  that 
provided  maximum  piece-wise  linearity  (as  determined  by  least  squares  methods) 
of  the  data  in  the  post-yield  stress  squared-strain  plane.  This  value  was 
0.822  GPa,  and  was  very  near  the  proportional  limit  stress.  Figure  10  shows 
the  data  intersect  the  strain  axis  at  0.395  percent,  (e  )  also  in  close  agree¬ 
ment  with  the  proportional  limit  value. 

An  important  feature  is  that  the  data  are  characterized  by  multiple 
slopes.  According  to  the  Bell  theory,  a  slope  change  can  occur  at  transition 
strains,  the  first  of  which,  from  Eqs.  (3)  and  (8)  with  N  =  18,  is  1.897  per¬ 
cent.  The  intersection  of  the  two  linear  segments  occurs  at  2.128  percent 
strain.  This  differs  from  prediction,  however,  scatter  in  the  values  of 
transition  strains  for  mild  steel  is  documented,*  an  exception  which  may  also 
apply  to  steel  alloys.  Note,  however,  that  no  change  in  slope  occurs  at  the 
second  transition  strain,  _  ^g* 

The  values  of  g,  the  deformation  moduli,  determined  from  the  least  squares 
fit  slopes  were  found  to  be  1.071  GPa  for  e  <  e  <  e^j  _  an<^  1*366  GPa  for 

e  >  e'  These  may  be  compared  with  the  values  predicted  by  Eq.  (6)  for 

r  =  3  (1.0868  GPa)  and  r  =  2  (1.3311  GPa)  mode  indices.  The  calculation  used 
the  experimentally  determined  values  of  E  and  v  and  assumed  the  melting  temper¬ 
ature  of  iron,  1808  degrees  Kelvin.21  The  experimentally  determined  values 


21  r.  c.  Weast ,  Ph.D,  ed. ,  Handbook  o£  Chemistri£  and  Physics,  56th  Ed. ,  The 
Chemical  Rubber  Co.,  Cleveland,  1975 ,  p  B-20. 
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Figure  9.  Log  wave  speed  vs.  log  dynamic  post-yield  strain  with  yield  strain  0.0046. 
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Figure  10. 


t-yield  stress  squared  vs.  strain  for 
5i-static  data;  yield  stress  0.822  GPa. 
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differ  by  less  than  two  percent  from  the  theoretical  values.  The  correlation 
of  theory  and  experiment  for  the  entire  strain  range  investigated  is  shown  in 
Figure  3,  where  the  parabola  coefficients  used  are  the  experimentally  deter¬ 
mined  values.  Table  3  summarizes  the  constitutive  equations  for  quasi-static 
loading . 

The  elementary  bar  theory  of  St.  Venant22  predicts  that  duration  of  impact 
should  be  independent  of  particle  velocity  and  dependent  only  upon  the  length 
of  the  bar  and  the  elastic  bar  wave  velocity.  The  St.  Venant  theory  is  com¬ 
pared  with  the  experimental  time  of  contact  data  in  Figure  5.  The  author 
believes  the  minimum  time  of  contact  and  the  bar  theory  should  coincide,  and 
the  discrepancy  is  due,  in  part,  to  experimental  technique.  In  a  comparison 
of  techniques  used  to  measure  time  of  specimen  contact  during  impact,  Filbey23 
shows  that  the  optical  technique  used  here  gives  values  approximately  seven 
percent  longer  than  those  obtained  by  electrical  contact  of  the  specimens  or 
piezocrystal  stress  measurements  at  the  impact  face.  Applying  this  correction 
to  all  the  data  points  brings  the  minimum  time  of  contact  in  fair  agreement 
with  the  bar  theory  prediction.  This  discrepancy  however,  has  no  effect  on 
the  particle  velocity  value  at  the  minimum  time  of  contact.  The  curves  in 
Figures  5  and  6  are  probable  curves  through  the  data  points  only,  and  no 
attempt  to  explain  their  particular  forms  has  been  made. 

The  dynamic  yield  stress,  found  to  be  16  percent  higher  than  the  quasi¬ 
static  value,  can  be  related  to  the  quasi-static  stress-strain  curve  by  the 
following  explanation.  The  stress  associated  with  the  first  transition  strain 
of  the  quasi-static  stress-strain  curve  is  0.954  GPa,  which  is  also  the  experi¬ 
mentally  determined  dynamic  yield  stress.  The  physical  significance  of  this 
phenomenon  is  that  a  material  accommodation  occurs  at  this  level  of  stress 
for  both  quasi-static  and  dynamic  loading.  It  is  interesting  that  the  quasi¬ 
static  deformation  modulus  increases  at  this  stress  while  the  dynamic  modulus 
decreases,  that  is,  yield  occurs. 

Elementary  elastic  wave  propagation  theory  predicts  that  elastic  strains 
propagate  with  the  elastic  bar  wave  velocity  CQ.  Studies  of  the  initial 
development  of  an  elastic  pulse  in  a  semi-infinite  aluminum  bar24  show  that 
the  theory  applies  at  large  distances  from  the  impact  end,  but  that  near  the 
impact  end,  the  arrival  time  of  low  level  elastic  strains  is  predicted  by  the 
dilatational  velocity.  The  experimental  data  of  Figure  8  show  this  to  be  true 
with  RHA  also.  The  maximum  axial  strain  which  may  propagate  with  the  dilata¬ 
tional  velocity  can  be  determined  by  calculating  the  dilatation  itself. 

Before  the  passage  of  the  shear  wave,  the  impacted  rod  experiences  only  a  wave 
of  dilatation  and  is  therefore  in  a  state  of  uniaxial  strain.  The  dilatation 
is  due  only  to  this  uniaxial  strain  until  the  shear  wave  arrives  and  produces 


22S.  Timoshenko  and  S.  N.  Goodier ,  Theory  of  Elasticity,  McGraw-Hill  Book  Co., 
Inc.,  New  York,  1951,  p  444. 

23 g.  L.  Filbey,  "Intense  Plastic  Waves,"  Ph.D.  dissertation.  The  Johns  Hopkins 
University ,  Baltimore,  Maryland,  1961. 

2V.  F.  Bell ,  "The  Initial  development  of  an  Elastic  Pulse  Propagating  in  a  _ 
Semi-Infinite  Bar,  "  BRL  Technical  Report  No.  6,  The  Johns  Hopkins  University, 
Baltimore,  Maryland,  November  1960. 
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the  Poisson  effect.  The  dilatation  A,  is  given  by 

A  =  P/K  (19) 

where  P  is  the  hydrostatic  pressure,  and  K  the  bulk  modulus.  The  hydrostatic 
pressure,  given  by  a^/3,  is  one-third  the  yield  stress  before  the  arrival  of 

the  slower,  higher  stress  plastic  waves.  The  computed  dilatation  is  0.20  per¬ 
cent  strain,  in  excellent  agreement  with  the  experimental  data  as  the  demarca¬ 
tion  between  the  two  wave  speeds  for  both  regions  of  the  rod,  one-to-two  and 
two-to-three  rod  diameters  from  the  impact  end.  Using  this  value  for  the 
dilatation  in  Eq.  (10)  a  dynamic  yield  stress  of  0.949  GPa  is  computed,  in 
close  agreement  with  that  computed  by  Eq.  (11) . 

The  progressive  decrease  of  wave  speeds  above  the  bar  wave  velocity  and 
simultaneous  increase  of  those  below  that  velocity  is  what  must  occur  for  the 
two-wave  structure  to  collapse  to  the  single  bar  wave  velocity.  An  estimate 
can  be  made  of  the  distance  from  the  impact  end  where  the  two  wave  structure 
converges  to  the  bar  wave.  Assume  that  the  initial  dilatational  and  shear 
velocities  have  accelerations  such  that  each  velocity  will  be  equal  to  the 
bar  wave  velocity  at  the  same  point.  For  this  to  occur 

(C1  -  ~  C2^  =  constant  (20) 

everywhere  along  the  rod,  where  C-^  and  are  the  mode  wave  speeds  of  strains 

below  and  above  the  dilatation  strain.  Using  the  values  of  and  to  com¬ 
pute  the  ratio  corresponding  to  the  impact  end,  the  constant  is  0.351.  The 
ratio  corresponding  to  a  position  two  and  one-half  rod  diameters  (midway 
between  positions  where’ arrival  times  were  measured)  is  0.358,  thus  establish¬ 
ing  the  validity  of  the  assumption  Eq.  (20).  Assuming  constant  accelerations, 
it  is  estimated  that  convergence  occurs  between  six  and  seven  diameter  lengths 
from  the  impact  end.  That  this  type  of  collapsing  two-wave  structure  could 
exist  is  shown  by  Bell25  for  annealed  aluminum  in  the  context  of  Truesdell ’ s25 
general  theory  of  waves  in  finite  elastic  strain. 

The  applicability  of  Bell’s  theory  to  the  quasi-static  deformation  data, 
and  the  experimental  evidence  that  wave  speeds  of  post-yield  strains  were 
governed  by  a  nearly  parabolic  stress-strain  relation  motivated  an  examination 
of  the  dynamic  deformation  data  in  the  context  of  the  theory.  The  exponent  a 
in  Eq.  (18)  was  assumed  equal  to  one-half,  and  the  wave  speed  expression  was 
linearized  with  respect  to  strain.  Wave  speed  to  the  minus  four  power  versus 
dynamic  post-yield  strain  was  found  to  be  a  linear  relationship,  as  Figure  11 
shows . 


25J.  F.  Bell ,  " Experiments  on  Large  Amplitude  Waves  in  Finite  Elastic  Strain , " 
Proceedings ,  IJJTAM  Symposium  on  Second  Order  Effects  in  Elasticity >  Plastic¬ 
ity  >  and  Fluid  Dynamics ,  Pergamon  Press ,  1964 ,  pp  173-186 . 

26C\  4.  Truesdelly  " General  and  Exact  Theory  of  Waves  in  Finite  Elastic  Strain , " 
Archy  Rational  Mech .  Anal.*  8j  No .  3>  1961 3  pp  263-352 . 
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Figure  11.  Linearized  wave  speed  vs.  dynamic  post-yield  strain 
data,  assuming  parabolic  response  function. 


The  value  of  the  deformation  modulus  3,  determined  from  the  slope  was 
1.092  1  .030  GPa.  This  is  within  two  percent  of  the  quasi-static  value  for 
this  strain  region,  and  very  close  to  theoretical  prediction.  It  cannot  be 
determined  whether  mode  index  changes  occur  at  the  transition  strains  for  this 
material  but  Bell’s  investigations  of  many  other  materials  show  that  mode  index 
changes  do  not  occur  in  dynamic  loading,  and  that  a  single  deformation  modulus 
3,  applies  as  long  as  stress  and  strain  are  monotonically  increasing. 4  Wave 
speed  versus  strain  is  plotted  in  Figure  12,  along  with  a  graph  of  the  theoret¬ 
ical  curve 

Cp  =  /e/2p  (e  -  1//4  (21) 


where  3  and  are  the  experimentally  determined  values. 

The  fact  that  plastic  wave  speeds  were  found  to  propagate  with  constant 
velocity  is  a  necessary,  but  not  sufficient  condition  for  the  applicability  of 
the  strain-rate  independent,  one-dimensional,  finite  amplitude  wave  propagation 
theory.  The  first  integral  relationship  [Eq.  (14)]  must  be  shown  to  apply. 
Particle  velocity  measurements  at  points  along  the  rod  were  not  made  during 
deformation,  so  only  the  terminal  condition  of  Eq.  (14) 


u 

max 


£max 

C(e) 


de 


(22) 


could  be  checked  with  experiment.  The  maximum  particle  velocity  can  be  parti¬ 
tioned  into  an  elastic  and  plastic  portion.  The  maximum  elastic  particle 
velocity  u^,  has  been  experimentally  determined,  and  Eq.  (21)  predicts  plastic 
wave  speeds,  thus: 


\ax  =  \  +  (4/3)  /g/2p 


(■emax  ey^ 


3/4 


(23) 


The  relatively  slow  velocity  of  the  plastic  wave  fronts  (compared  to 
elastic  waves) ,  coupled  with  finite  specimen  lengths  restricted  the  propaga¬ 
tion  distance  of  large  strains  because  their  amplitude  was  reduced  by  the 
reflected  elastic  precursor.  However,  low  level  plastic  strains  did  propagate 
larger  distances  due  to  their  higher  velocity,  and  a  plateau  of  maximum  strain 
was  observed  on  rods  struck  at  low  velocity.  This  can  be  seen  in  Figure  7, 
Test  20,  where  the  maximum  strain  at  50  and  75  mm  from  the  impact  end  is  the 
same.  This  maximum  strain  is  0.515  percent  and  is  predicted  by  Eq.  (23)  for 
the  maximum  particle  velocity  for  this  test.  The  maximum  strain  for  the 
higher  impact  velocities  did  propagate  unreduced  at  least  the  first  rod  diam¬ 
eter  length  from  the  impact  end.  It  was  assumed,  as  a  first  approximation, 
that  the  magnitude  of  the  unloading  strain  had  been  equal  to  yield  strain,  and 
that  the  maximum  strain  was  the  sum  of  the  yield  strain  and  the  permanent,  or 
residual  strain. 

e  =  e  +  e  (24) 

max  y  w 


That  the  residual  strain  measured  at  the  first  diameter  position  corresponded 
closely  with  the  prediction  of  Eqs.  (23)  and  (24)  is  shown  in  Figure  13  for 
rods  with  a  maximum  particle  velocity  of  0.040  mm/ps.  The  predicted  residual 
strain  is  1.678  percent,  and  the  size  of  the  symbol  indicates  the  scatter  of 
the  individual  measurements. 
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Figure  12.  Wave  speed  vs.  compressive  strain  compared  with 
prediction  of  parabolic  response  function. 
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Figure  13.  Residual  strain  vs.  distance  from  impact 
end  for  rods  struck  at  80  m/s. 


Strain  maxima  at  one  rod  diameter  length  from  the  impact  end  were  computed 
by  Eq.  (24),  and  are  plotted  in  Figure  14  (circles)  against  the  maximum  particle 
velocity  of  the  respective  test.  A  graph  of  Eq.  (23)  with  as  the  indepen¬ 

dent  variable  (solid  line)  shows  the  correlation  of  the  theory  and  experiment. 
The  experimentally  determined  value  of  3  (1.092  GPa)  was  used,  and  this  single 
value  of  deformation  modulus  predicts  the  strain  maxima  even  beyond  the  first 
transition  strain  of  1.960  percent.  This  is  consistent  with  the  observation 
that  mode  index  changes  do  not  occur  in  dynamic  loading.4 

The  apparent  departure  of  theory  from  experiment  in  predicting  maximum 
strains  for  the  higher  impact  velocities  is  due  to  the  progressively  decreasing 
validity  of  Eq.  (24)  .  Using  momentum  considerations,  Lee2'  shows  that  the 
reflected  elastic  precursor  is  gradually  absorbed  as  it  encounters  the  propa¬ 
gating  plastic  wave  front.  This  is  mathematically  described  by 

6a  =  -  (pCQ  u^  -  a^)/ 2  (25) 


where  the  subscript  b  refers  to  values  before  the  arrival  of  the  unloading 
wave,  and  6a  is  the  change  in  stress  at  the  plastic  wave  front  after  passage 
of  the  unloading  wave.  Assuming  linear  elastic  unloading,  the  unloading  strain 
is  given  by 


e  i 

uni 


6a/E 


(26) 


As  particle  velocity  and  stress  increase,  the  magnitude  of  unloading  strain 
approaches  zero.  If  it  is  assumed  that  u^  is  given  by  Eq.  (23)  and  a^  by 

Eq.  (27)  below  for  the  same  level  of  maximum  strain,  the  strain  maxima  are 
given  by  the  squares  in  Figure  14,  which  lie  in  closer  agreement  to  theoretical 
prediction . 

This  correlation  of  theory  and  experiment  in  predicting  maximum  strains 
from  wave  speeds  is  deemed  sufficient  that  the  strain-rate  independent,  one¬ 
dimensional,  finite  amplitude  wave  propagation  theory  applies.  Thus,  dynamic 
stress  is  given  by  Eq.  (15),  which  becomes 

a  =  ay  +  3(e  -  Ey)1/2  (27) 

when  the  wave  speed  expression  [Eq.  (21)]  is  the  integrand,  and  the  subscript 
y  refers  to  dynamic  yield  values.  Furthermore,  the  dynamic  deformation  modulus 
3  is  given  by  Eq.  (6)  for  a  mode  index  r  =  3.  The  quasi-static  and  dynamic 
stress-strain  curves  are  compared  in  Figure  15. 


VI.  SUMMARY  AND  CONCLUSIONS 

Experimentally  determined  constitutive  equations,  which  describe  the  quasi¬ 
static  and  dynamic  response  of  thirty-eight  millimeter  thick  RHA  plate  to  uni¬ 
axial  compression,  are  presented.  The  quasi-static  post-yield  stress-strain 


27  E.  H.  Lee ,  "A  Boundary  Value  Problem  in  the  Theory  of  Plastic  Wave  Propaga¬ 
tion , "  Quart.  Appl.  Math ,  10,  No.  4,  1953,  pp  335-346. 
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Figure  14.  Maximum  strain  vs.  maximum  particle  velocity  compared  with  finite  amplitude 
wave  propagation  theory  using  parabolic  response  function. 
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Figure  15.  Comparison  of  static  and  dynamic  stress- 
strain  curves  for  38  mm  RHA. 
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response  function,  in  nominal  form,  is  parabolic  with  the  origin  at  the  yield 
point,  corresponding  to  the  intermediate  region  of  Bell's  general  theory  of 
parabolic  plasticity.  Two  discrete  deformation  moduli  characterize  the 
response  function.  The  change  in  modulus  occurs  at  the  first  transition 
strain  of  the  general  theory,  provided  that  the  finite  yield  strain  is 
accounted  for.  The  moduli  themselves  compare  favorably  with  those  given  by  a 
single,  quantized,  temperature  dependent  relationship  based  on  the  value  of 
the  isotropic  elastic  shear  modulus  at  absolute  zero  temperature. 

The  dynamic  yield  stress  is  higher  than  that  measured  in  quasi-static 
loading,  but  is  related  to  the  latter  in  that  the  dynamic  yield  stress  and 
the  quasi-static  stress  at  the  first  transition  strain  are  equal.  This 
suggests  a  material  accommodation  that  is  stress,  not  strain  dependent.  The 
dynamic  stress-strain  response  function  was  deduced  by  using  a  one-dimensional, 
strain-rate  independent,  finite  amplitude  wave  propagation  theory,  the  applica¬ 
bility  of  which  was  shown  experimentally.  The  dynamic  response  function  also 
is  parabolic,  with  its  origin  at  the  dynamic  yield  point.  A  single  deformation 
modulus,  nearly  equal  to  the  modulus  of  the  initial  portion  of  the  quasi-static 
stress-strain  curve,  is  applicable.  This  further  relates  the  quasi-static  and 
dynamic  response  functions.  Wave  propagation  velocities  and  strain  maxima  near 
the  impact  end  are  well  predicted  by  the  wave  propagation  theory  when  the 
parabolic  response  function  is  inserted.  However,  strain  arrival  times  and 
propagation  distance  of  maximum  strain  are  not,  indicating  that  the  maximum 
stress  is  not  initiated  as  a  step  function  at  the  impact  face. 

Elastic  waves  near  the  impact  end  do  not  propagate  with  a  single  velocity, 
but  with  a  two-wave  structure.  Strains  below  the  maximum  expected  dilatation 
propagate  with  the  dilatational  velocity,  while  elastic  strains  greater  than 
the  maximum  dilatation  propagate  with  the  shear  velocity.  The  trend  of  the 
data  indicate  that  this  two-wave  structure  degenerates  to  that  of  the  bar  wave 
at  distances  far  from  the  impact  end. 

Although  the  strain  range  investigated  is  limited  to  six  percent,  deforma¬ 
tion  of  RHA  by  uniaxial  compression  is  explained  by  a  general  theory  of  plas¬ 
ticity  for  metal  alloys.4  The  functional  form  of  the  response  function  is 
preserved  over  a  strain  rate  range  of  seven  orders  of  magnitude,  and  the 
increase  in  dynamic  stress  is  attributable  to  the  elastic,  not  plastic  region. 
Additional  experiments  to  extend  and  bridge  the  strain  and  strain-rate  ranges 
investigated  would  provide  the  data  needed  to  extend  the  applicability  of  the 
Bell  theory  to  stresses  and  impact  velocities  of  ballistic  interest. 
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LIST  OF  SYMBOLS 


a  exponent,  dimensionless 

e  coefficient  of  restitution,  dimensionless 

r  deformation  modulus  mode  index,  dimensionless 

t  time,  ys 

u  displacement,  mm 

B  deformation  modulus  constant,  0.0280 

o 

C  wave  velocity,  mm/ys 

C  elastic  bar  wave  velocity,  mm/ys 

o 

C1  elastic  dilatational  wave  velocity,  mm/ys  • 

C2  elastic  shear  wave  velocity,  mm/ys 

D  specimen  diameter,  mm 

E  Young's  modulus,  GPa 

H  power  function  coefficient,  GPa 

K  bulk  modulus,  GPa 

L  specimen  length,  mm 

N  integer  constant,  dimensionless 

P  hydrostatic  pressure,  GPa 

R  strain  gage  resistance,  ohms 

R  initial  resistance  of  unstrained  strain  gage,  ohms 

o 

5  deformation  modulus,  GPa 

T  ambient  material  temperature,  degrees  Kelvin 

a 

Tm  material  melting  temperature,  degrees  Kelvin 

X  Lagrangian  coordinate,  mm 

3  plastic  deformation  modulus,  GPa 

6  change  in  quantity  operator 

e  nominal  engineering  strain,  mm/mm 
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'N 


'N 


'uni 


X 

VI 

V 

p 

0 
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plastic  transition  strain,  mm/mm 

total  transition  strain,  mm/mm 

unloading  strain,  mm/mm 

final,  permanent  strain,  mm/mm 

Lame's  constant,  GPa 

shear  modulus,  GPa 

Poisson's  ratio,  dimensionless 

mass  density,  g/cm^ 

nominal  engineering  stress,  GPa 
3  3 

dilatation,  mm  /mm 


Subscripts 

b  value  of  quantity  before  unloading 

max  maximum  value  of  quantity 

p  quantity  in  plastic  condition 

y  value  of  quantity  at  yield 


Superscripts 

first  derivative  with  respect  to  time 
_  statistical  mode  value  of  quantity 
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